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SUMMARY 


Some  of  the  trends  in  the  design  of  wings  and  control  surfaces 
for  high  speed  vehicles  have  brought  about  a  need  for  modifications  in 
the  methods  of  analysis  for  stress  distribution  and  flexibility.  A 
few  characteristics  typical  of  such  structures  are  low  aspect  ratio, 
large  angle  of  sweep,  thinness,  and  the  use  of  thick  skins,  sometimes 
stiffened  by  closely  spaced  integral  spars  and  sometimes  forming  sand¬ 
wich  structure. 

The  present  report  outlines  some  modifications  that  can  be  applied 
to  the  matrix  force  method  when  the  effects  of  some  of  these  character¬ 
istics  are  important.  These  refinements  do  not,  of  course,  affect  the 
application  of  the  method  to  more  conventional  structures.  It  has 
already  proven  to  be  adequate  without  the  refinements  for  most  of  the 
structures  handled  to  date. 

The  report  is  concerned  with  applications  of  the  matrix  force 
method  to  wing  structures  idealized  as  consisting  of  capstrips  carry¬ 
ing  axial  load  only,  interconnected  by  shear  panels  which  can  carry* 
shear  only.  The  refinements  being  introduced  are  for  the  purpose  of 
taking  into  account  the  Poisson'1  s  ratio  effect  in  the  case  of  thick 
skin  wings  and  for  improving  upon  the  handling  of  the  stress  coupling 
effect  in  the  case  of  non- orthogonal  ribs  and  spars.  A  system  of  * 

choosing  unit  redundant  load  distributions  is  also  discussed.  These 
refinements  have  all  been  treated  in  some  form  or  another  by  other 
authors.  However,  previous  methods  of  accounting  for  the  Poisson's 
ratio  effect  apparently  have,  except  in  certain  special  cases,  neglected 
the  fact  that  while  the  strain  is  the  same  in  a  spar  cap,  say,  and  in  the 
adjacent  cover  material  in  the  same  direction,  the  stresses  may  be  dif¬ 
ferent  because  of  a  state  of  bi -axial  stress  in  the  cover.  A  way  of 
taking  this  into  account  has  been  included.  The  question  of  whether  this 
difference  can  be  significant  has  not  yet  been  resolved,  and  it  may  turn 
out  to  be  negligible  In  all  practical  cases. 

These  modifications  are  reviewed  in  theory,  and  the  formulas  and 
procedures  required  for  practical  application  are  developed  and  pre¬ 
sented  in  a  table.  The  Poisson's  ratio  refinement  is  then  applied 
to  the  case  of  a  homogeneous  laterally  loaded  flat  plate,  and  excellent 
agreement  with  an  available  plate  theory  solution  is  shown. 
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SYMBOLS 


— aiii 

^  X.’  Oy*  ^Xjr 

E,  G 

V 


X 

b,J 


U 


a ,  h 
Mx'  1 

A 


w 

(j  V,(f  S,  f-TB 


=  symbols  defined  in  ref.  1 

=  stress  components  referred  to  rectangular  coordinates 
=  Young’s  Mod-ulus,  Shear  Modulus 
=  Poisson's  Ratio 

=  ratios  defined  where  used  in  the  text 
=  skin  thickness 

=  stress  components  referred  to  oblique  coordinates 
=  sweep  angle 

=  cell  dimensions  in  rib  and  spar  directions  respectively 
=  strain  energy 

=  side  and  depth  respectively  of  square  plate 

=  bending  moments  per  unit  length  of  sections  of  plate 
perpendicular  to  the  x  -  and  y  -  axes  respectively 

=  lateral  displacement  of  plate 

=  intensity  of  a  distributed  load 

=  stress  components  referred  to  rectangular  coordinates 
with  an  axis  in  the  spar  direction. 
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DISCUSSION  OF  THEORY 


Introduction 

The  report  discusses  first,  for  clarity,  the  Poisson’s  ratio  refine¬ 
ment  applied  to  an  unswept  wing.  Next,  a  swept  wing  made  up  of  identi¬ 
cal,  parallelogram  cells  is  treated  and  finally  the  more  practical  case 
of  a  swept  wing  which  has  taper  and  non-parallel  ribs  and  spars.  The 
formulas  for  the  necessary  flexibility  influence  coefficients  for  this, 
the  most  general  case,  are  given  in  a  table.  Formulas  for  the  simpler 
cases  can  be  deduced  from  these.  Although  the  derivations  of  the 
formulas  may  appear  complicated,  the  formulas  themselves  and  the  instruc¬ 
tions  for  using  them  are  summarized  compactly  in  the  table. 


Idealization  -  Unswept  Case 

Consider  the  problem  of  analyzing  an  unswept,  thick  skin,  multi-cell 
box  beam  as  shown  in  Figure  la.  Following  the  conventional  matrix  force 
method  approach,  one  lumps  the  covers  into  equivalent  capstrips*  over 
the  spars  and  ribs;  the  contributions  of  the  ribs  and  spars  themselves 
are  calculated  and  added  in.  If  appropriate,  intermediate  shear  lag 
members  are  introduced  as  well.  Connecting  the  various  spar  and  rib 
caps  are  panels  as sumed  to  be  capable  of  carrying  shear  loads  only.  Thus, 
the  axial  loads  in  the  capstrips  vary  linearly  from  one  panel  point  to 
the  next. 

Based  on  this  idealization,  the  required  load  distributions  in  the 
statically  determinate  structure  can  be  obtained  in  any  of  a  number  of 
ways.  A  recommended  method  is  discussed  later.  The  flexibility  matrix 
is  calculated  from  the  formulas  to  be  given  and  the  matrix  operations 
are  then  carried  out  in  the  usual  manner. 

For  the  derivation  of  the  member  flexibilities  for  a  structure  of 
this  type,  it  is  convenient  to  augment  the  lumped  bar  and  shear  panel 
idealization  with  a  further  assumption.  It  is  assumed,  for  the  purpose 
of  accurately  expressing  the  strain  energy  in  terms  of  the  member  loads 
that  the  stiffening  effect  of  the  discrete  ribs  and  spars  would  not  be 
appreciably  altered  if  they  were  spread  out  over  the  wing  covers;  in 
other  words,  that  the  stiffened  cover  can  be  adequately  represented  by 
an  orthotropic  plate  with  properly  assigned  deflectional  characteristics. 
This  assumption  is  employed  in  the  following  manner. 

Consider  a  typical  rectangular  segment  of  the  rib-  and  spar- stiffened 
cover,  as  shown  in  Figure  lb.  The  strain  energy  in  such  a  segment  is 
simply  the  sum  of  the  skin  and  stiffener  energies.  If  the  strain  in 
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DISCUSSION  OF  THEORY 


Idealization  -  Unswept  Case  (Continued) 

each  direction  (spar  and  rib)  is  assumed  not  to  vary  over  the  segment, 
the  energy  can  be  written  as 

U  =  h[^s2  +C*r2  +  1  Tsr2]i  4  t  +  aspar  4  ^spa/ 

ft  ^  2  (1) 

+  arib  z"  ^rib  , 

where  (j  ,  (j  and  T  are  the  skin  stresses,  (see  Fig.  lc),  and 


- -  '  -  ax 

and  t  are  “the  dimensions  of  the  segment . 

It  has  been  taken  into  account  in  this  expression  that  in  a  stiffened 
cover  the  stress  in  the  cap  is  different  from  the  axial  stress  in  the 
skin  in  the  cap  direction.  This  is  because  there  is  a  Poisson's  ratio 
effect  in  the  skin  but  not  in  the  cap. 

To  simplify  the  strain  energy  expression,  the  symbols  rs  and  rr  are 
introduced  and  defined  as  the  ratios  of  stiffener  area  to  total  cross 
sectional  area  of  the  segment  in  the  spar  and  xjb  directions  respectively. 
Also  defined  are  the  average  stresses  0S  and  Q  r,  as  the  total  load 
divided  by  the  cross  sectional  area  in  each  of  the  two  directions.  Sub¬ 
stituting  these  quantities  into  (l),  the  strain  energy  expression  can 
be  manipulated  into  the  form 


2E  [  es  l-re  vs  rr  l-rr  ■ 

+  I  rsr2]^s  4  t 


r.C 
r  sr 


2V<Je  (Jr 


where 


]J  rr^  -  rs) 
°ss  =  1  '  1  -y*  r  r 


V 2  rs(l  -  rr ) 


n  =  i  -  1 - -sx - 1 

rr  1  -y  rsrr 

,  .2 

c  - 1  +JZ-E6S: - 

sr  1  1  -L/2  rsrs 
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Idealization  -  Unswept  Case  (Continued) 

In  this  expression  the' C’s  account  for  the  difference  in  stress 
between  the  skin  and  stiffeners.  If  this  difference  is  neglected,  each 
of  these  quantities  reduces  to  unity.  As  previously  stated,  the  stif¬ 
feners  are  assumed  to  be  spaced  closely  enough  to  permit  treating  the 
integrally  stiffened  cover  as  an  orthotrppic  plate.  With  this  assump¬ 
tion  made,  the  strain  energy  of  an  infinitesimal  element  of  the  equiv¬ 
alent  orthotropic  plate  is  simply  expression  (2)  with  XQ  replaced  by 
dJ?s  and  J?r  replaced  by  dflr.  It  is  now  possible  to  evaluate  the  member 

flexibilities  by  integrating  this  expression  over  the  entire  structure. 

« 

JIhe  integration  is  carried  out  as  follows.  Referring  to  Figure  2, 
the  (j  *  and  q”s2  terms  are  integrated  piece-wise  over  volumes  of  material 
that  are  the  same  as  were  lumped  into  bars  panels  previously.  For 
example,  in  Fig.  2a  the  integration  of  the  0r2  term  is  carried_out  over 
the  cross-hatched  area.  Fig.  2b  shows  the  area  over  which  the  (j  s2  term 
is  integrated.  With  the  stress  varying  linearly  between  panel  points, 
the  resulting  flexibility  influence  coefficients  are  the  same  as  those 
obtained  in  ref.  1  except  for  the  "C"  factors  which  take  into  account 
the  difference  in  stress  between  skin  and  stiffener.  The  7*sr  term 
integrated  over  a  panel  as  shown  in  2c  gives  the  same  expression  for 
shear  flexibility  as  ref.  1.  As  suggested  in  ref.  3,  the  assumption 
is  made  that  since  the  contribution  of  the  Poisson's  ratio  term.  Os Ov> 
to  the  total  strain  energy  is  somewhat  smaller  than  that  of  the  direct 
stress  terms,  a  less  accurate  approximation  is  sufficient  for  it. 
Accordingly,  in  approximating  the  integral  of  this  term  over  the  struc¬ 
ture,  it  is  assumed  that  the  stresses  (jQ  aud  (jr  are  constant  over  areas 
of  which  the  cross-hatched  portion  of  Fig.  2d  is  typical,  and  are 
to  be  the  stress  values  at  the  panel  point  corner.  An  example  would  be 
member  loads  qg  and  qo  in  Fig.  2e.  The  resulting  flexibility  influence 
coefficient  formulas  for  one  cover  become 


a  1,1  -  2C<1,2  -a2,2  ■  Crr  • 

^3,3  =  2(^3,^  =  Css 

^  =  (l/Et)bi!  (2  +  2 1))  or  A/Gt 


(3) 


CX  o  o  ~  -  C  (l/Et)LR  R 
'-a2,3  sr'  '  ^  s  r 

where  the  subscripts  appended  to  the  CX  * s  refer  to  the  member  loads  shown 

in  Fig,.  2e.  R  is  the  ratio  of  skin  area  (cross-hatched  side  only,' Fig. 2e) 
£ 
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DISCUSSION  OF  THEORY 


Idealization  -  Unswept  Case  (Continued) 

to  total  luuiped  spar  area  in  a  section  normal  to  the  spar,  and  similarly 
p  -the  ratio  of  skin  area  to  total  lumped  rib  area  in  a  section  normal 
to  the  rib.  The  average  skin  thickness  of  the  area  of  integration  in 
each  case  is  t.  By  employing  flexibility  influence  coefficients  of  this 
type  with  the  conventional  approach,  the  Poisson’s  ratio  effect  can  be 
adequately  taken  into  ac count. 

It  can  be  shown  that  in  the  present  case,  where  there  is  no  sweep, 
the  "C"  factors  will  usually  be  within  a  few  percent  of  unity  and  in 
no  event  more  than  10$  different  from  unity.  Since  the  Poisson  s  ratio 
effect  is,  itself,  responsible  for  only  a  small  variation  in  the  stress 
distribution,  it  can  be  concluded  that  consistency  will  be  maintained, 
if  these  terms  are  taken  equal  to  unity.  This  assumption  and  supporting 
argument  are  also  made  in  ref .  8  where  a  similar  problem  has  been  treated. 

It  may  be  noted  that  the  member  flexibility  formulas  given  in  ref.  2 
for  swept  wings  reduce  to  the  formulas  (2)  if  there  is  no  sweep  and  the 
effect  of  the  C's  is  neglected. 

Idealization  *»  Swept  Case 

As  the  next  step  in  the  development  of  the  formulas  required  for 
the  analysis  of  a  practical  swept  wing  structure,  consider  the  case  of 
a  swept  wing  with  uniform  depth,  spars  uniformly  swept  back  and  identical 
ribs  all  parallel  to  the  vehicle  centerline.  (See  Fig.  4a)  The  infini¬ 
tesimal  element  of  interest  in  this  case  is  a  parallelogram,  and  it  is 
convenient  to  use  oblique  stresses  as  shown  in  Figure  3«  These  are  the 
same  as  were  employed  by  Hemp  in  ref.  7*  (The  formulas  derived  in  this 
section  have  been  shown  to  be  in  agreement  with  the  work  of  ref.  7*  > 

Again  defining  average  stresses 


3*.- 


aB  sec  ^ 


and  d  - - 

wr  By.  sec^ 


the  strain  energy  expression  for  the  element  can  be  written.  This  expres¬ 
sion  and  subsequent  algebraic  manipulation  leading  to  the  flexibility 
formulas  is  given  in  Appendix  A.  The  energy  expression  differs  from  l 2 ) 
in  two  important  respects.  First _there  is  the  introduction  of  the  terns 
which  couple  the  direct  stresses  (jB  end  Qr  with  the  shear  stress  TSr* 

As  pointed  out  in  ref.  3,  these  are  the  terms  which  account  for  the 
peaking  of  the  spanwise  stress  in  the  covers  of  swept  wings  in  the  vlcln- 
ity  of  the  rear  beam  at  the  wing  root.  (ref.  1  allows  for  the  7 sr_ps 
portion  of  the  interaction;  the  TBr(jr  *>art  is  omiivfced;  however.  ef .  1 
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Idealization  -  Swept  Case  (Continued) 

7-  2 

sr  contribution,  but  leaves  out  the  Poisson's 
ratio  term  Os  Or  altogether. )  Secondly,  there  is  the  appearance  of 
the  parameter  fj^  (  =Ucob2  A“  sin2  \  )  which  has  replaced  Poisson's 
ratio,  1/ ,  in  eqn.  (2).  As  a  result,  unless. the  sweep  angle  is  small, 
the  "C"  factors  may  become  significantly  different  from  unity  and  may 
have  to  be  retained  in  the  formulas. 

The  recommended  procedure  is  to  again  lump  the  covers  into  equiva¬ 
lent  capstrips  over  the  ribs  and  spars  with  intermediate  shear  lag 
members  as  required.  By  integrating  the  strain  energy  expression  for 
the  infinitesimal  oblique  element  over  the  structure,  making  similar 
assumptions  as  in  the  unswept  case,  the  member  flexibilities  can  again 
be  determined.  The  Q  2  and  ljr 2  terms  yield  the  same  flexibilities 
as  before.  The  Tsr2  term  is  now  multiplied  by  a  function  of  the  sweep 
angle . 

In  the  case  of  the  interaction  terms,  it  is  again  sufficiently 
accurate  to  assume  that  £fg  and  0r  are  constant  over  areas  of  which 
the  shaded  portion  of  Fig.  4a  is  representative,  and  that  they  are 
equal  to  the  values  at  the  corner  member  loads,  as  q£  and  qo  in  Fig.  4b. 

The  flexibility  influence  coefficient  formulas  for  the' swept  wing  with 
identical  parallelogram  cells  are  therefore  (referring  to  Fig.  4b): 

an  =  2Ct12=a22  =  crr  -b/3alE 


&  33  =  2  “  ^44  =  Css  ’  ^/3a3E 

a  55  =  C  *(l/Bt)bJe  cos  X  (2  +  2  4  tan2  A) 

a  23  =  ~Csr  '  (1/Et)RsRr  ( V -  tan2  A)  cos  A 


&25  =  tan  A 

OC  35  =  CTs  *  (VEt)R/  tan  A 
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DISCUSSION  OF  THEORY 


Idealization  -  Swept  Case  (Continued) 

In  the  preceding  Rg  is  again  the  ratio  of  skin  area  (cross-hatched 
side  only,  Fig.  4b)  to  total  lumped  spar  area  in  a  section  normal  to 
the  spar,  and  R^.  is  again  similarly  defined.  An  outline  of  the  deriva¬ 
tion  of  these  formulas  is  given  in  Appendix  A.  (These  formulas  were 
given,  without  detailed  discussion  and  with  the  "C"  factors  neglected, 
in  ref.  2. ) 


Idealization  -  Non-Parallel  Spars  and  Non-Parallel  Ribs 

The  preceding  two  sections  have  discussed  cases  where  the  wing 
exhibited  uniformity  of  stiffness  and  regularity  of  planform  that  would 
rarely  be  found  in  actual  practice.  The  more  practical  formulas  and 
procedure  given  in  this  section,  however,  can  best  be  developed  with 
the  preceding  formulas  as  a  basis.  In  most  actual  cases,  the  spars 
of  a  wing  structure  are  not  parallel;  if  in  addition  the  wing  is  swept, 
some  of  the  ribs  are  frequently  in  the  free  stream  direction,  while 
others  are  normal  to  the  structural  axis.  In  such  cases,  the  5s2  and 
Qr2  strain  energy  terms  are  still  adequately  accounted  for  by  the 
lumped  stringer  approach.  However,  the  stringers  will  generally  be 
tapered  and  the  flexibility  expressions  will  consequently  be  preceded 
by  the  already  tabulated  factors,  "4"  (See  ref.  l),  which  take  the  taper 
into  account.  The  panels  may  turn  out  to  be  of  any  general  quadrilateral 
shapes.  For  these,  the  strain  energy  term  can  be  handled  by  the 

expression  for  quadrilateral  shear  panels  derived  as  an  extension  of 
the  work  of  Garvey  (ref.  6),  and  given  in  Appendix  B.  This  leaves  only 
the  Poisson’s  ratio  and  sweep  coupling  terms  unaccounted  for. 

It  is  suggested  that  these  terms  be  handled  as  follows.  Fig.  5 
shows  a  general  quadrilateral  plate,  bounded  by  spars  and  ribs.  The 
cross-hatching  indicate®  the  area  for  which  the  coupling  terms  are  to 
be  approximated.  Because  of  the  non-parallelogram  shape,  the  shear 
flows  kj  qc;  and  kg  q^  are  unequal.  Imagine  that  the  given  panel  is 
replaced  by  the  parallelogram  panel  which  is  indicated  by  the  dashed 
lines,  and  take  the  shear  flow  acting  on  it  to  be  kqc ,  where  k  =  ■§-  (k-l+k2). 
The  Poisson's  ratio  terms  may  now  be  approximated  bythe  fourth  of 
eqns.(ij-).  The  sveep  coupling  terms  may  be  approximated  by  the  last  2 
of  eqns.  (4),  modified  for  the  k  factor. 

A  suggested  scheme  for  calculating  averaged  values  of  the  C's  is 
given  inthe  table  on  the  following  page,  as  are  the  flexibility  influence 
coefficient  formulas  for  a  swept  wing  of  arbitrary  planform  and  taper. 
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Idealization  -  Non-Parallel  Spars  and  Non-Parallel  Ribs  (Continued) 

The  treatment  of  the  coupling  terns  is  admittedly  a  crude  approxi¬ 
mation.  However,  these  terms  exert  a  decidedly  lesser  influence  upon 
structural  deflections  and  stress  distributions  than  do  the  direct 
stress  terms,  and  thus  the  rougher  approximation  is  actually  consistent 
with  the  rest  of  the  analysis.  The  trend  is  indicated  in  the  example 
given  later  on.  In  the  case  of  various  practical  structures  that  have 
been  analyzed,  a  comparison  between  the  results  obtained  with  these 
coupling  terms  included  and  those  with  the  terms  omitted  altogether, 
further  substantiates  this. 

After  the  redundant  internal  load  distribution  has  been  determined 
using  the  above  method,  the  skin  axial  and  shear  stresses  in  rectangular 
coordinates,  and  the  spar  and  rib  stresses,  can  be  calculated  by  the 
stress  transformation  formulas  given  in  Appendix  C. 


Selection  of  Redundants  and  Applied  load  Distributions 

A  matrix  formulation  practically  the  same  as  in  ref .  1  has  been  used 
by  Grzedzielski  (ref.  3)  and  by  Argris  (ref.  4)  in  conjunction  with  a 
scheme  which  greatly  simplifies  and  systematizes  the  determination  of 
the  number  of  redundants  and  the  choice  of  redundant  distributions  in 
multi-cell  structures.  A  discussion  of  the  theory  need  not  be  repeated 
in  this  report,  but  the  results  are  essentially  that: 

1.  The  number  of  redundancies  (excluding  shear-lag  members)  in  a 
multi-cell  structure,  with  no  discontinuities  and  not  including 
support  redundancies,  is  equal  to  the  number  of  internal  walls, 
and 

** 

2.  A  simple  and  systematic  set  of  redundant  distributions  (again 
excluding  shear  lag  members)  can  be  formed  by  applying  equal 
opposite  warping  groups  to  each  pair  of  adjacent  cells  in  the 
structure  at  their  common  wall  (see  Fig.  6). 

If  shear  lag  members  are  present,  additional  locally  self-balancing 
redundant  distributions  can  be  easily  devised  to  handle  them,  (see,  for 
example,  pg.  10,  ref.  4) 

As  pointed  out  in  ref.  4,  the  applied  load  distributions, 

(statically  determinate  applied  load  distributions  by  the  customary  term¬ 
inology),  may  be  selected  from  any  distributions  that  equilibrate  the 
applied  load.  They  should,  therefore,  be  made  as  simple  as  possible. 

(f  or  example :  cantilever  beams,  simply  supported  beams,  distributions 
according  to  engineering  beam  theory,  etc.)  Strictly  speaking,  these 
distributions  need  not  necessarily  be  statically  determinate  but  are 
usually  selected  as  such  for  convenience. 
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EXAMPLE 


Description 

A  simple  example  Illustrating  the  results  of  the  use  of  some  of  the 
discussed  techniques  is  presented  in  this  section.  The  structure  is  a 
Square,  simply  supported  solid  plate  with  side  a  and  depth  h,  subjected 
to  uniform  lateral  load.  This  particular  example  was  chosen  for  two 
reasons;  an  exact  solution  is  available  for  comparison,  ref.  5>  and  it 
represents  the  limiting  case  of  a  thick  skinned,  high  density  core  wing. 

If  the  recommended  semi-monocoque  analysis  approach  works  here,  it  should 
be  even  better  for  the  more  usual  wing  structure. 

The  plate  is  idealized  as  a  2  x  2  grid,  shown  . schematic ally  in  Fig.  7a, 
and  also  as  a  4  x  4  and  a  6  x  6  grid.  In  addition  to  solutions  based 
upon  the  recommendations  of  the  present  report,  each  case  is  also  solved 
omitting  the  Poisson's  ratio  term.  In  this  way  the  relative  importance 
of  this  effect  can  be  qualitatively  evaluated. 

In  determining  the  equivalent  capstrip  areas,  the  familiar  technique 
of  maintaining  identical  bending  stiffness  is  observed;  see  Fig.  7b. 

The  cover  plates  are  determined  in  the  same  manner;  thus  the  equivalent 
thickness  is  taken  as  h/6.  The  vertical  shear  webs  are  assumed  to  be 
infinitely  stiff;  this  is  in  agreement  with  plate  theory  assumptions. 

,  The  applied  loads  used  to  simulate  the  uniformly  distributed  loading 
are  concentrated  forces  acting  at  the  interior  grid  points.  In  order  to 
hold  down  the  complexity  of  the  calculations,  only  unit  loads  exhibiting 
double  symmetry  are  considered.  Thus  it  becomes  necessary  to  analyze  only 
one  eighth  of  the  plate,  giving  it,  of  course,  the  proper  boundary  condi¬ 
tions.  The  numbering  system  for  designating  the  applied  loads  is  given 
in  Figure  Jc. 


Results 


A  comparison  of  bending  moments  and  deflections  at  the  center  of  the 
plate  is  given  in  Table  1  for  each  case  investigated.  The  exact  solution 
of  ref.  5  is  used  as  the  basis.  Bending  moment  and  deflection  distribu¬ 
tions  along  the  symmetry  axes  are  also  presented  for  the  exact  solution 
in  Fig.  8,  together  with  the  corresponding  semi-monocoque  solution  results. 
Table  2  gives  the  calculated  flexibility  influence  coefficients. 

An  examination  of  Table  2  indicates  that  the  central  deflection  due 
to  a  unit  central  load  agrees  very  well  with  the  exact  result,  the  differ¬ 
ence  varying  from  15$  for  the  2x2  grid,  to  2$  for  the  6x6. 
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EXAMPLE 


Results  (Continued) 

The  central  deflection  due  to  uniformly  distributed  load.  Table  1, 
shows  a  more  extreme  variation,  the  difference  ranging  from  kOjo  for  the 
2  x  2  to  h$>  for  the  6x6.  It  is  felt  that  the  bcrfo  value  is  the  result 
of  the  excessively  crude  representation  of  the  distributed  load  by  a 
single  central  force.  On  the  other  hand,  the  central  bending  moment 
for  the  2x2  appears  to  be  artificially  close  to  the  exact  value*  As 
the  number  of  grid  points  increases,  the  representation  of  the  distributed 
load  by  concentrated  loads  should  become  better  and  better,  hnd  both 
central  deflection  and  central  moment  agreement  reflect  this. 

Tables  1  and  2  also  Indicate  that  the  influence  of  the  Poisson's 
ratio  effect  is  quite  important  •  in  this  example,  being  approximately 
twice  as  much  so  in  the  case  of  deflections  as  stresses.  The  effect 
should  diminish  somewhat  in  the  case  of  wing  structures  which 
have  concentrations  of  capstrip  material  over  the  webs  of  the  spars 
and  ribs. 

In  conclusion,  it  is  felt  that  the  example  demonstrates  that  the 
method  of  accounting  for  the  Poisson's  ratio  effect  as  recommended  in 
this  report  is  quite  satisfactory.  Furthermore,  it  would  appear  .reason¬ 
able  to  inffer  that  the  recommended  method  for  handling  the  somewhat  similar 
sweep-coupling  effect  should  work  equally  well. 
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APPENDIX  A 

Strain  Energy  in  a  Stiffened  Swept  Plate 


v 


\  \  Rib  \ 

- v  n - 

\  \  Spar 
\  \ 


VT*  C/y 


V  <4 


V  =  Poisson’s  ratio 

A  =  sweep  angle 

fJL  ~  L/cos2  A  -  sin2  A 

qg,  =  mem^er  loa(3-s  spar  and  rib  directions,  lb. 

q  =  shear  flow,  lh/in 

t  =  plate  thickness,  in. 

as  _  =  skin  area  plus  spar  area  in  section  normal  to  spar 

ay  =  skin  area  plus  rib  area  in  section  normal  to  rib 

E,  G  =  Young  *  s  Modulus,  Shear  Modulus  lb/ in^ 

R_  =  ratio  of  skin  area  to  skin  area  plus  spar  area  in  section 

8  normal  to  spar 

Ry  =  ratio  of  skin  area  to  skin  area  plus  spar  area  in  section 

normal  to  rib 

rg  =  ratio  of  spar  area  to  skin  area  plus  spar  area  in  section 

normal  to  spar  ’ 

rr  =  ratio  of  rib  area  to  skin  area  plus  rib  area  in  section 

normal  to  rib 


Css,Crr>Csr 


=  constants  in  strain  energy  expression 


W°rr 
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APPENDIX  A 


Discussion 


The  uniform  thickness  swept  plate  shown  in  the  illustration  is 
reinforced  by  an  integral  spar  and  rib.  The  stresses  are  applied  such 
that  the  axial  strain  in  the  spar  direction  is  uniform  throughout  plate ^ 
and  spar,  and  likewise  in  the  rib  direction.  A  shear  stress  7  is  applied 
to  the  plate  as  well.  Because  of  the  Poisson's  ratio  effect,  the  spar 
stress  is  different  from  the  corresponding  stress  in  the  plate,  and 
similarly  for  the  rib. 

Using  the  standard  stress-strain  relations,  an  expression  for  the 
strain  energy  in  the  stiffened  plate  is  obtained,  as  follows. 

Strain  Energy  =  Plate  Energy  +  Spar  Energy  +  Rib  Energy 

Strain  Energy  =  jtbA  co6  A  (O^2  +  O^^sec2  A-  2fJi  (j s  (jT  sec  A 

+  kT  (  (3  )  tan  X  sec  A  +  T  2(if  +  ^  tan2  A  ) 

sr  s  r  sr 

+  a  i  (  (j  )2  +  a  b((l r.?v)2')  (l) 

spar  ^  '“'spar  rib  rib  j 


Defining  average  stresses 


V s  =  ag  sec^  ’  -  ar  sec^ 


and  q/t 


and  using  the  definitions  of  q£J  qr  and  q,  the  following  are  derived. ' 

G"s  -  B  [<7  +Wrs  S’r  ‘  2  sln  Ars  tea  *  l)rJ 
<Jr  *  B  ^rr  ^  -  2  sl”  Arr  ^rs  +  l)T] 


T«r  =  r 
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APPENDIX  A 


Discussion  (Continued) 
and 


^spar  =  [(1  -ffrr^B  ‘  M  ^  “  rs}^r  +  2  sin^(l  '  rs}^rr  +  l)f 

C^rib  =  ['A^1  '  +  (1  “M2rs^r  +  2  sin^(l  "  ^(^s  +l)7V 


where 


spar 

IW  =  r-r - t-*- - 

s  bt  cos  1  +  a 


r  =  rib 

r  ~  KT^jTT 


D  =  1  -//  r  r 
M  s  r 

Substituting  expressions  (3)  and  (4)  into  (l),  the  strain  energy 
for  the  segment  shown  in  Fig.  4a  can  be  written  as 


.  -1 —  ft 
1  -  r  ur 


sec 

Ac  -A- 

-  O'2  +  c 

[  SB  1  -  ] 

r  s  rr 

s 

sr 

•  2Mcrscfr  + 

_  ,E  2 

Vr*  (g  cos 

*  k  sin  A  5" 

+  C  *  4  sin 

sr 

s 

rT 

2  )  \«-2 


where 


C  =l-i 
ss 


C  =  1  -/ 
rr 


C  -  1  + 
sr 


rr(l  -  r£ 


q  -  i  .S/rAy.a^r  -L£fl 

rr  d 


c  =  1  + 
sr 


C  =  1  + 

rr 


,rs(l  +  jLL  ri 


-rAl  +  JLcrf 
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APPENDIX  B 


Flexibility  Influence  Coefficient  for  a  Quadrilateral  Shear  Panel 
From  Reference  0  ~ 


1.  General  Quadrilateral  Panel 


The  flexibility  influence  coefficient  for  the  general  quadrilateral 
panel  shown  above  is,  from  ref.  6: 

rsf  -  v  T  1  G  2F  7 

^all  Gt  *  E  TaTc’7(b-d)(a+b+c +d)J 

where 

F  =  (a+b)  +  ^(a^+b^)  +  y(a^+b^)  log  (a-fb) 

+  (c+d)  +  -|(c3+d3)  +  |.(c5+d5)  log  (c+d) 

-  (b+c)  +  |-(b3+c3)  +  -^(b5+c5)  log  (b+c) 

-  (d+a)  +  j.(d3+a3)  +  f(d5+a5)  log  (d+a) 

+  ^(a2-c2)(b3-d3)  +  £(b2-d2)(a3-c3)  -  |(a-c )(b4-d^) 

*  y(b-d)(a  -c  )  -  J-(a-c)(b-d)(a+b+c+d) 
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APPEHDIX  C 

Calculation  of  Skin,  Spar  and  Rib  Stresses 


The  use  of  oblique  stress  components  was  recommended  to  simplify  the 
calculation  of  member  flexibilities  in  the  case  of  swept  wings.  For  the 
purpose  of  designing  the  wing,  the  maximum  stresses  are  required.  These 
can  be  calculated  from  the  idealized  structure  member  loads  by  the  pro¬ 
cedure  to  follow. 

To  get  the  skin  stress  components  referred  to  a  rectangular  system 
of  coordinates  as  shown  in  Fig.  4c,  the  following  formulas  (derived  from 
expressions  (3)  and  (4)  in  Appendix  A  and  the  formulas  for  transformation 
of  stress  components  from  oblique  to  rectangular  coordinates)  are  applicable. 

»  mm 

(j'B  1  +/isin2A  /it  sin2 A  -2  ,sln2A)+l 

=  1 

(j'T  1"/ik  (1  cos2  A  cos2A  -2((X+l)  sin  A  cosA 

T  sin  A  cos  A  sinAcosA  -2(JJ+l)  sin2  A  +  1  q/t 

In  the  above  equation  the  difference  in  stress  between  skin  and  stif¬ 
feners  (due  to  Poisson’s  ratio)  has  been  neglected.  Also,  in  evaluating 
the  stresses  at  a  point  on' the  wing  skin  it  will  be  necessary  to  average 
the  neighboring  values  of  qs/ag,  qr/ a^.  and  q/t  in  some  appropriate  manner. 

To  obtain  the  magnitude  and  direction  of  the  maximum  skin  stresses, 
it  is  necessary  only  to  apply  to  the  primed  stress  components  the  stan¬ 
dard  formulas  for  principal  stresses  in  rectangular  coordinates. 

4 

The  spar  and  rib  stresses,  tinder  the  same  assumption,  are  given  by 

C\  =  q  /a  and  (j  =  q/a 
spar  us'  s  '■'rib  mr  r 
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TABLE  2 

FLEXIBILITY  INFLUENCE  COEFFICIENTS 

iy.i  r 

.  ■  ,  .1265 


.  ,  .1265*** 
1  .1(073** 

I  .1651* 

2x2  Grid 
Idealization 


.1221 

.0788 

.0523 

.1603 

.1010 

.067 9 

2  ‘ 

.0567 

.0391*1 

.071*0 

.05050 

3 

.03053 

.01*09 

1 

.1265 

.121*6 

.1581* 

1  -  •  ■  1 

C  -  .  ’• 

.1007 

.1267 

.0833 

*1053 

{  "i 

f>  ■  ■ ,  _ 

.0532 

.0671 

.01*51 

.057,2 

.0251 

.0321 

: 

2 

.  .  ;  ..  .  .  '.r  .  ■ 

r  ■ 

.0861 

.1090 

.0729 

.0919 

.01*77 

.0603 

.0l*0l* 

.0511 

.0226 

.0286 

3 

- 

.o6i*o 

.0811 

.01*16 

.0526 

.036I* 

.01*60 

.0208 

.0263 

. . 

jh  • 

.0301* 

.0389 

.0252 

.0317 

.011*1 

1  .0175 

5 

[’  •;  ■ 

.0228 

.0290 

.0139 

.0171* 

6 

• 

.0108 

.011*1 

6  x  6  Grid  Idealization 

i*“  *5®  *bOTe  "“‘‘rices  is  the  displacement  (in  non- 
EK/f  )  •*  a  P°tat  ■  in  the  plate  due  to  a  nnit  load  at  a 
point  n*  The  designation  of  points  is  given  in  Figure  7c. 

***  Plate  theory  solution.  Ref.  5 

. **  Redundant  analysis,  Poisson *s  ratio  term  included 

*  Redundant  analysis,  Poisson* s  ratio  term  omitted 
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